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Within the Dirac theory of the electronic properties of graphene, smoothly varying lattice strain 
affects the Dirac carriers through a synthetic gauge field. For static lattice strain, the gauge field 
induces a synthetic magnetic field which is known to suppress weak localization corrections by a 
dynamical breaking of time reversal symmetry. When the lattice strain is time dependent, as in con- 
nection with phononic excitations, the gauge field becomes time dependent and the synthetic vector 
potential is also associated with an electric field. In this paper, we show that this synthetic elec- 
tric field has observable consequences. We find that the Joule heating associated with the currents 
driven by the synthetic electric field dominates the intrinsic damping, caused by the electron-phonon 
interaction, of many acoustic phonon modes of graphene and metallic carbon nanotubes when in- 
cluding the effects of disorder and Coulomb interactions. Several important consequences follow 
from the observation that by time-reversal symmetry, the synthetic electric field associated with 
the vector potential has opposite signs for the two valleys. First, this implies that the synthetic 
electric field drives charge-neutral valley currents and is therefore unaffected by screening. This 
frequently makes the effects of the synthetic vector potential more relevant than a competing effect 
of the scalar deformation potential which has a much larger bare coupling constant. Second, valley 
currents decay by electron-electron scattering (valley Coulomb drag) which causes interesting tem- 
perature dependence of the damping rates. While our theory pertains first and foremost to metallic 
systems such as doped graphene and metallic carbon nanotubes, the underlying mechanisms should 
also be relevant for semiconducting carbon nanotubes when they are doped. 



I. INTRODUCTION 

It is one of the remarkable aspects of the Dirac descrip- 
tion of the low-energy electronic properties of graphene^'- 
that both disorder and lattice strain give rise to synthetic 
gauge fieldsi^i^i^i^ii If the synthetic vector potential A(r) 
is due to static disorder or static lattice distortions, its 
presence affects the electronic dynamics through an effec- 
tive magnetic field eB = V x A which, by time reversal 
symmetry, points in opposite directions at the two Dirac 
points of the electronic dispersion. Consequences of this 
effective magnetic field have been widely studied in the 
literaturei^i^i^ d°di Additional physics arises when the 
synthetic vector potential is caused by time- dependent 
distortions such as phonons. In this case, A becomes 
time dependent and will generate not only a magnetic 
but also an effective electric field eE = —dA/dt. 

In this paper, we show that these synthetic electric 
fields have observable consequences. Consider a low- 
energy phonon mode of graphene or carbon nanotubes 
(CNT). The phonon is associated with a synthetic elec- 
tric field which, when the system is metallic, drives cur- 
rents. The dissipation (Joule heating) associated with 
these currents causes damping of the phonon mode. We 
find that frequently, the synthetic electric fields are di- 
rectly responsible for the damping of phonon modes in 
metallic carbon nanotubes and graphene. In the clean 
limit, this damping mechanism is equivalent to dissipa- 
tion by electron- hole pair creation. In fact, we find that 
we reproduce corresponding recent results for the radial 
breathing mode of clean carbon nanotubes, Approach- 
ing the problem from the point of view of the synthetic 
electric fields allows us to calculate damping rates includ- 



ing the effects of disorder and electron-electron interac- 
tions which we find to be significant. Moreover, we find 
appreciable damping rates even for those phonon modes 
of carbon nanotubes for which damping by electron-hole 
pair creation is not effective due to the large discrepancy 
between the electron and the phonon velocities. As a re- 
sult, we expect damping by the electron-phonon interac- 
tion to dominate over other mechanisms such as phonon- 
phonon couplingi^ over a wide range of parameters. 

Searching for observable consequences of the synthetic 
electric field is complicated by the fact that lattice dis- 
tortions do not only induce a vector potential A but also 
a scalar potential^ (f>{r, t) which leads to an additional 
electric field —Vcj). Within a tight-binding approach, the 
synthetic vector potential is associated with changes in 
the hopping amplitude caused by changes in the bond 
length. The scalar potential is a deformation potential 
which arises from local dilation or compression of the lat- 
tice. In fact, estimates'^ suggest that the bare coupling 
constant of the scalar potential is about an order of mag- 
nitude larger than that of the vector potential. 

Nevertheless, we find that the damping of most (but 
not all) low-energy phonon modes of graphene and car- 
bon nanotubes is dominated by the vector potential. This 
is a consequence of two important qualitative differences 
between the electric fields associated with the scalar and 
vector potentials: (i) While the scalar field affects both 
valleys in the same manner, the sign of the vector po- 
tential is opposite for the two valleys. As a result, the 
currents driven by the electric fields are charge currents 
for the scalar potential, but valley currents for the vector 
potential. Therefore, the valley electric fields due to the 
vector potential will not be subject to screening as they 
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do not induce any charge densities. Since in contrast, the 
electric fields arising from the scalar deformation poten- 
tial are screened, this may significantly affect the relative 
importance of the two electric fields when the system is 
metallic, (ii) The scalar potential is necessarily associ- 
ated with a longitudinal electric field whose direction is 
parallel to the wavevector. In contrast, the vector poten- 
tial generates a synthetic electric field which in general 
has both longitudinal and transverse components. We 
find that over a wide range of wavevectors, the transverse 
conductivity (and hence the corresponding dissipation) is 
significantly larger than the longitudinal conductivity. 

The fact that the vector potential drives valley currents 
has another important consequence. While by momen- 
tum conservation, electron-electron interactions cannot 
induce a decay of charge currents, they do lead to a de- 
cay of valley currents by a process which is analogous 
to spin Coulomb dragJ^ii^ii^ We find that this valley 
Coulomb drag mechanism leads to interesting tempera- 
ture dependence of the phonon damping which in many 
cases is decreasing as temperature increases. 

The damping mechanisms which we discuss apply most 
directly to systems with metallic behavior, i.e., metallic 
carbon nanotubes and graphene. For this reason, we will 
not explicitly discuss semiconducting CNT in this paper. 
It has recently been suggestedi^ that the intrinsic vibra- 
tional damping in semiconducting CNT is dominated by 
nonlinear elastic effects. However, it should be kept in 
mind that in practice, many semiconducting CNT are 
doped and can hence exhibit a finite conductivity. In 
this case, the damping mechanisms discussed in this pa- 
per may well be relevant as well. 

It is found experimentally that the vibrational re- 
laxation time of the radial breathing mode is remark- 
ably long, of the order of several nanoseconds ?il This 
follows from scanning-tunneling-microscope-based trans- 
port measurements, in which the tunneling current ex- 
cites the radial breathing mode, which enables vibron- 
absorption processes at temperatures far below the 
phonon frequency. We find that our results for the vi- 
brational relaxation time of the radial breathing mode 
are consistent with these experiments. 

Identifying parameters for which long vibrational re- 
laxation times can be realized in CNT is also of much 
interest from the perspective of nanoelectromechanical 
systems made of carbon nanotube s^^i^''''^^i^^'^°i^^'^^ and 
graphenei^"^'^^ For example, if vibrational relaxation 
times are sufficiently long, these systems provide access 
to a regime in which the transport current drives the vi- 
brational mode far out of thermal equilibrium. For strong 
electron-vibron coupling, this nonequilibrium regime can 
be characterized by self-similar avalanche transport i^!^ 
One of the most promising systems to observe this ef- 
fect are suspended carbon nanotube quantum dots for 
which strong electron-vibron coupling and the associated 
Franck-Condon blockade have recently been confirmed 
experimentally.^^ 

This paper is organized as follows. In Sec. |TT1 we in- 



troduce the strain-induced scalar and vector potentials 
fSec. lII A)) and derive a general expression for the relation 
between the phonon damping rate and the conductivity 
(Sec. IIIB[) . We also illustrate the basic physical picture 
in the context of the radial breathing modes of carbon 
nanotubes (Sec. Ill C)) . Our approach allows us to include 
the effects of disorder (Sec. IIIII) and of electron-electron 
interactions (Sec. IIV[) on phonon damping. Finally we 
conclude in Sec. |Vl Some calculational details are rele- 
gated to an appendix. 



II. BASIC FORMALISM 

A. Strain-induced vector and scalar potentials 

Elastic strains couple to electrons in graphene and 
carbon nanotubes by inducing effective scalar and 
vector potentials into the low-energy electronic Dirac 
Hamiltonian,"^ 



n 



£v^>\v) {VF [i(T,{d, - Ai)] + 0} l'(r) 



(1) 



where ^'(r) is the spinor which defines the electron, 
is a Pauli matrix, vp is the Fermi velocity, and we show 
the Hamiltonian defined for one of the two, K and K' , 
valleys in the Brillouin zone of graphene. Here we chose 
coordinates such that the unit vectors of the graphene 
lattice can be written as 



ai = aV3 I —^x + -jTiy 
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where a « 1.4A denotes the bond length. Both the vector 
and the scalar potential can be expressed in terms of the 
strain tensor Uy . The scalar potential is determined by 
the trace of the strain tensor, (j){v,t) — gniu^^ + Uy 
where gn is estimated to be of order 20-30eV in Ref. 
The form of the vector potential is essentially fixed by 
symmetry to be^i^ 



A(r) = Mf 

2a V '^xx 



xy 



(3) 



where /3 = d\og{t)/d\og{a) w 2 — 3. Strictly speaking, 
this expression has a small uncertainty in the prefactor 
since the theory of elasticity may not accurately describe 
the displacements within the unit cell. 

In nanotubes, the natural coordinate system, defined 
by the nanotube axis, is rotated by an angle 9 with re- 
spect to the coordinate axes of graphene defined above. 
Choosing the CNT axis as the x-axis and the direction 
around the tube as the y-axis, the angle 9 is given by 
cot 9 ~ V3 "7'" for (n,m) carbon nanotubes. It takes the 
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value = (6* = 7r/2) for zigzag (armchair) CNT. The 
vector potential takes the form 



A(r) 



h(3 

2a 



D{3d) 



yy 



(4) 



where D{39) is a rotation matrix. Note that in this equa- 
tion, also the strain tensor is given in the rotated coor- 
dinate system. 

The acoustic phonon modes of carbon nanotubes and 
graphene can be described within standard elasticity the- 
ory. For graphene, the elastic Lagrangian density for the 
strain tensor Uij and the out-of-plane displacement h(r, t) 
takes the form 



stretch 



end 



(5) 



with 



T 
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Stretch — M^ij ^ 



(6) 



in terms of the 2D mass density poi the Lame coefficients 
/i and A characterizing the in-plane rigidity of the lattice, 
and the bending rigidity k. The same Lagrangian also 
applies to carbon nanotubes of radius i?, when replacing 
the bending energy by 
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The strain tensor takes the form Uij = {l/2)[diUj+djUi + 
{dih){djh)] for graphene and 

ox 
du 
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for carbon nanotubes. Here, u denotes the displacements 
within the graphene sheet and h the displacement in the 
perpendicular (for CNT: radial) direction. 



B. Phonon damping 

As argued in the Introduction, there is a close rela- 
tion between phonon damping and the conductivity ten- 
sor. This relation can be obtained formally by computing 
the shift Alo in the phonon frequency due the electron- 
phonon coupling. The damping rate is then given by 
r = 2ImAci;. Consider first the case in which the phonon 



is associated with a vector potential A{r,t). We can ex- 
press the vector potential in terms of phonon creation 
and annihilation operators, bq and 6^, respectively, 



A(r,t) = ^^ylqexp(iqr) 
q 

X (^bqexp{-iu!qt) + b^_qexp{iu!qt)^ (9) 



with 
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Here, we defined the matrix 
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and a unit vector p describing the mode polarization 
in Ux^Uy, and /i-direction, which takes the form p = 
{l,0,iXqxR/{2^ + A)) for the longitudinal stretching 
mode, p = (0, 1,0) for the twist mode, and p — (0,0, 1) 
for the radial breathing mode [as obtained from the 
Euler-Lagrange equations for the elastic Lagrangian Eq. 
([S])]. The surface area of the nanotube is denoted by 
51 — {2i:R)L and the mode dispersion by ojq. The cor- 
responding results for the acoustic modes of graphene 
follow by taking R^ oo and setting 6 = 0. (The flexural 
modes of of graphene will be discussed separately below.) 

In second order perturbation theory, we then find for 
the damping rate of a phonon with wavevector q. 
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Using the Kubo formula, we can express Fq in terms 
of the dissipative (real and symmetric) contribution 
crs(q, cj) to the conductivity tensor. 
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Here, v = vpcr denotes the velocity operator. Note that 
we employ the two-dimensional conductivity unless ex- 
plicitly stated otherwise. Thus, we can express the damp- 
ing rate as 

^''=2^^-(^'-)(-^- + ^^^)^-^- + i^'>)- 

(14) 

Here, we have also included the effect of the scalar po- 
tential d> and defined 
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in terms of the vector 

M^ = itq,,iqy,l/R). 



(16) 



Note that we have included the effects of screening in the 
contribution originating from the scalar potential, (f (q) 
denotes the Coulomb interaction and n(q, us) the polar- 
ization operator.) 



C. Basic physical picture: Damping of the radial 
breathing mode in armchair CNT 

Before we embark on a systematic investigation of 
phonon damping based on the expressions derived above, 
we would like to illustrate the basic physics in the context 
of the radial breathing mode of armchair carbon nan- 
otubes. We will do this in the context of a semiclassical 
approach which clearly brings out the physics and com- 
plements the more quantum mechanical approach taken 
in the remainder of the paper. This example also shows 
that our approach yields damping rates which are of the 
order of those observed in experiment. 

For the radial breathing mode of carbon nanotubcs, 
only the radial displacement h is nonzero (i.e., u = 0) in 
the long- wavelength limit. According to Eq. the cor- 
responding strain tensor takes the form Uxx = u^y = 
while Uyy = h/R. Using that for armchair carbon nan- 
otubes 9 = 7r/2, the synthetic vector potential takes the 
form 



A(r,i) 



h{r,t)/R\ 
2a [ J 
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Similarly, we find for the scalar deformation potential 
<j>{r,t) = guh{r,t)/ R. In the long wavelength limit 
q ^ 0, the electric field — V(/) associated with the scalar 
potential vanishes. At the same time, the valley electric 
field eE = —dA/dt originating from the vector potential 
remains finite because the RBM dispersion tends to a fi- 
nite frequency ujb = [(2/^ + A)//Oo^^]^^^ for q ^ 0. Thus, 
the synthetic vector potential gives the dominant con- 
tribution to the damping of the radial breathing mode. 
Note also that the electric field is pointing in the direction 
along the nanotube axis. 

For clean armchair carbon nanotubes, the damping of 
the radial breathing mode vanishes to leading order,.!^ 
since the contribution of the synthetic vector poten- 
tial to the Hamiltonian is proportional to ax and hence 
commutes with the unperturbed Hamiltonian. However, 
since armchair carbon nanotubes are metallic, there will 
be damping in the presence of disorder. Describing the 
optical conductivity of carbon nanotubes within a Drude 
model, we have 



<Jdc 



1 
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in terms of the mean free time r. From the Einstein re- 
lation, the (one-dimensional) dc conductivity takes the 



form adc = where £ = vpT denotes the electronic 

mean free path and N counts the spin and valley degener- 
acy. We can now compute the damping rate of the radial 
breathing mode, F, by equating the time derivative of the 
elastic energy £ with the Joule heating associated with 
Recr(ti;). 

The elastic energy of the radial breathing mode (per 
unit length) can be readily obtained from the elastic La- 
grangian so that we find 

f ^ I [{2nR)pou;%\h{.B)n 



— Recr(aj_B) 



(19) 



Inserting e|E(w)| - ^ 
field, we find a damping rate F of 



^^\h{Lu) \ for the synthetic electric 



8TT'^{poa^)R^l + {ujBT) 
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(20) 



An interesting aspect of this expression is that the damp- 
ing rate falls off with the third power of the nanotube 
radius. This may be useful guidance to enter into 
the regime of current-driven nonequilibrium in nanoelec- 
tromechanical devices. 

Numerical estimates of Eq. yield relaxation times 
of the order of nanoseconds for an elastic mean free 
path of 1 fim and a diameter of order Inm. This is of 
the same magnitude as the relaxation times observed in 
experiment 



III. DISORDER EFFECTS ON PHONON 
DAMPING 

A. Carbon nanotubes 

Due to the weakness of screening in one dimension, the 
scalar potential remains relevant in carbon nanotubes. 
Indeed, one readily estimates that in the absence of 
screening, the ratio of the electric fields due to vector 
and scalar potential is of the order 
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(21) 



for phonon modes with a linear dispersion u = cq. At 
the same time, the suppression of the scalar potential by 
screening 

[1 + i;(q)H(q, c.)]-! [1 + — Hl/qR)]-^ (22) 

involves, for realistic values of q, a factor smaller than 
but still of order one. Here, v denotes the electronic den- 
sity of states and R the radius of the nanotube. As a 
result, we conclude that the damping of the longitudinal 
stretching mode is dominated by the effects of the scalar 
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potential. In contrast, the damping of the radial breath- 
ing mode (whose dispersion approaches a finite frequency 
as g — > 0) will be dominated by the vector potential. In 
the following, we will only discuss the dominant dissi- 
pation channel for the various phonon modes of carbon 
nanotubes. An evaluation of the subdominant dissipa- 
tion channel would, however, be straight-forward. 



1. Radial breathing mode 




qo 



i/i 



logq 



In order to make contact with the literature,— it is in- 
structive to start with the damping of the radial breath- 
ing mode in clean carbon nanotubes. One readily estab- 
lishes that only the j/-component of the vector potential 
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(23) 



contributes to dissipation. Evaluating Eq. p3p . we ob- 
tain for the corresponding component of the dissipative 
conductivity tensor 



Ne?VF 



(24) 



Inserting Eqs. (p3)) and (I24p into the general expression 
Eq. p4|l , we recover for the damping rate of the RBM of 
clean carbon nanotubesii 



RBM 



8{poa'^)R^Lu 



cos^{W)e{u) -2vFkF)- (25) 



Specifically, this expression confirms the absence of 
damping for armchair CNT where 6 — 7r/2. This shows 
that our approach is equivalent to standard approaches 
to damping of phonon modes. 

The strength of making the relation with the synthetic 
electric fields lies in allowing us to go beyond the limit 
of clean and non-interacting samples. Including disorder 
for the armchair CNT, there will also be damping due 
to electric fields pointing along the CNT. Again using 
a Drude expression for the frequency-dependent conduc- 
tivity, one recovers from Eq. the semiclassical result 
given above in Eq. (PO)) . 



2. Longitudinal stretching mode 

At long wavelengths, the longitudinal stretching mode 
of carbon nanotubes has a much lower frequency lo = 
cq than the radial breathing mode and we can make a 
diffusive ansatz for the dynamic conductivity, 



cr(q,tj) 



-lujadc 



(26) 



which is valid when q£ 1 and D = vf^/^ is the dif- 
fusion coefficient. Note that the diffusion pole in this 



FIG. 1: (Color online). Sketch of the damping rate as function 
of wavevector for the longitudinal stretching mode of carbon 
nanotubes. 



expression is protected by charge conservation since the 
scalar potential drives ordinary charge currents. The dis- 
sipative conductivity is then given by 



Re(T(q, uj) 



{Dq^y 



(27) 



The precise q dependence of the conductivity depends 
sensitively on the relative magnitudes of w = cq and 
Dq^. To clearly bring out the q dependence of the 
damping rate, we consider the long-wavelength regime 
q ^ qoj dominated by w, and the short- wavelength 
regime q ^ qo, dominated by Dq^, separately. Here, 
the characteristic wavevector dividing between these two 
regimes is given by 



„ c 1 
go -2 — -. 
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Evaluating the scalar potential Eq. (|15p for the longitu- 
dinal stretching mode and inserting the resulting expres- 
sion into Eq. (|14[) . we obtain the damping rate 
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Here, we defined the renormalized coupling constant 
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which, strictly speaking, still includes a weak logarithmic 
q dependence. 

Eq. (|29p predicts that the damping rate increases 
quadratically with q for long wavelengths and saturates 
to a q-independent constant for shorter wavelengths q S> 
go- This behavior is sketched in Fig. [T] It is interest- 
ing to compare Fq to the mode frequency Wq in order to 
see whether the longitudinal stretching mode can become 
overdamped for some region of wavevectors. Clearly, the 
ratio Fq/wq is maximal for q = qq. Remarkably, the 
maximal value of this ratio becomes independent of the 
elastic mean free path and thus quite universal. 
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FIG. 2: (Color online) Orientation of the synthetic electric 
field at the two K-points of the dispersion (red arrows) with 
respect to the phonon wavevector q for (a) the purely trans- 
verse and (b) the purely longitudinal situation. As shown 
in Eqs. (|33p - (I36|) . the synthetic electric field is in general 
neither purely parallel nor transverse. 



electric field is purely longitudinal, the weaker bare cou- 
pling of the vector potential is partially offset by the fact 
that dissipation by longitudinal valley currents is very 
sensitive to (and strongly enhanced by) disorder-induced 
intervalley scattering. (It is interesting to note that in 
the context of multi-valley semiconductors, the impor- 
tance of intervalley scattering in acoustic attenuation is 
actually known for many decades. ^^^) 

Using Eq. we can readily obtain the longitudinal 
and transverse components of the vector potential. Con- 
sistent with the lattice symmetry, we find 



hp 2hn . . 

2a y po^q 
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for the longitudinal phonon 
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(34) 

(35) 
(36) 



Inserting numbers typical of carbon nanotubes into this 
expression, one finds values of order 0.1 showing that 
even though the longitudinal stretching mode remains 
underdamped at all wavelengths, damping can be quite 
significant. 



B. Graphene 

1. Longitudinal and transverse acoustic (in-plane) phonons 

The analysis of the previous section can be extended to 
phonon modes of graphene. However, phonon damping 
in (doped) graphene will typically be dominated by the 
vector potential. At long wavelengths, this is a conse- 
quence of the much stronger screening in two dimensions 
than in the carbon nanotube setting, 
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l + i;(q)n(q,c^) 



qtf 



(32) 



where qxF = 2Tre^v denotes the Thomas-Fermi wavevec- 
tor of graphene. (For a density of approximately 
lO^^cm^^, the screening length is of the order of tens 
of nm.) At short phonon wavelengths (relative to the 
elastic mean free path), this emerges from the fact that 
unlike —Vcj), E = —dA/dt will almost always have a sig- 
nificant transverse component (relative to the wavevec- 
tor q, cp. Fig. [2]).— It turns out that in this range of 
wavevectors, the transverse conductivity (and hence the 
associated damping) is much larger than the longitudi- 
nal conductivity. In addition, even when the synthetic 



for the transverse phonons. Here, 6q denotes the angle 
between the direction of the wavevector q and the a;-axis. 

At finite q, the dissipative conductivity becomes a 
symmetric tensor which is diagonal in a coordinate sys- 
tem whose axes are parallel and perpendicular to the 
wavevector q. 



a{q,uj) 



CT||(q,w) 
crj_{q,uj) 



(37) 



For doped graphene in the diffusive regime ql <^ 1, 
the longitudinal and transverse conductivities can be ob- 
tained from hydrodynamic equations for charge densities 
and currents which follow from a Boltzmann equation. 
One finds the continuity equations 



dni 
~dt 
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as well as Ohm's laws 
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(41) 



Here, the indices 1 and 2 label the valleys, v denotes 
the density of states, and 1/r and 1/ry are the intra- 
valley and inter-valley scattering rates due to disorder, 
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FIG. 3: (Color online). Schematic dependence of the damp- 
ing rate Fq on wavevector q of acoustic phonons in graphene 
when the synthetic electric field is purely longitudinal, cp. 
Eqs. (|50p and (|5H) . This applies when the wavevector points 
in a narrow cone around the zigag (armchair) direction for the 
longitudinal (transverse) acoustic phonon. Left: qo > ly^. 
Right: go < 



respectively. We left out the effects of induced electric 
fields from these equations because screening does not 
affect the valley odd channel. We can now obtain the 
valley-odd conductivity relating E_ = Ei — E2 to j_ — 
ji ~j2 by taking the difference between the two equations 
()40|) and ()41|) . combined with the continuity equation in 
the valley-odd channel. 

For the longitudinal conductivity, we find 



cr||(q,w) = 



(42) 



The reason for the cutoff of the diffusion pole by the inter- 
valley scattering rate 1 /ry is that there is no conservation 
law associated with valley currents which is analogous to 
charge conservation for ordinary charge currents. Despite 
the weakness of intervalley scattering (originating from 
atomic-scale defects) compared to intravalley scattering, 
it is actually rather important to include l/ry. From Eq. 
(|42| . the dissipative conductivity takes the form 



ReCT||(q,tj) = 
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TV 
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TV 



)]o-dc 



(43) 



While it is evidently possible to work with this complete 
expression, it is more instructive to analyze the various 
limiting cases. Except for the scale qo = {2c/vf){1/£) 
introduced above, this expression involves the interval- 
ley scattering length iy = [^DtvY^"^ as a second length 
scale. When q£v ^ 1, one finds that ReCT||(q, w) ~ crdc- 
In the opposite limit q£v 3> 1, one obtains 

f o-dc g < go ^ ^ 

Re CT„ (q, uj) ~ i i+{qoiv? (44) 

Transverse fields do not induce any charge densities and 
therefore, we find 

CT_L(q,t^) = o-dc (45) 
for the transverse conductivity at any q Xjl. 
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FIG. 4: (Color online). Schematic dependence of the damp- 
ing rate Fq on wavevector q for acoustic phonons in graphene 
when the tranverse synthetic electric field is dominant, cp. 
Eq. (|49|l . This applies for all directions of q, except when 
the wavevector points in a narrow cone around the zigag 
(armchair) direction for the longitudinal (transverse) acoustic 
phonon. 



For doped graphene in the ballistic regime ^ 1, we 
can obtain the dynamic conductivity (T(q, lo) from the 
Boltzmann equation. Heuristically, we can obtain fT(q, w) 
(up to numerical prefactors) by making the replacements 
Dq^ ^ vpq and l/£ ^ q in the diffusive results Eqs. ([^5)1 
and (|45)l . A more formal derivation is relegated to the 
Appendix. For transverse electric fields, disorder- induced 
intervalley scattering is irrelevant (since no valley charge 
densities are induced) so that the transverse conductivity 
becomes 



CT_L(q,t^) = 



(46) 



for uj <^ vpq- It is interesting to note that this trans- 
verse conductivity cr^ also controls the attenuation of 
surface acoustic waves in the fractional quantum Hall ef- 
fect at Landau level filling factor v = 1/2.^^ In contrast, 
intervalley scattering is important for the longitudinal 
conductivity when q^tv ^ 1; 



while for q^ly 3> 1, we find 
cr||(q,w) = 



1^ 1^^777 



e vuj 



vpq 



(47) 



(48) 



for all g ^ l/£. One readily establishes that a-\\/a± ~ 
[c/vfY' ^ 10^'* for sufficiently large 9, confirming that 
the damping will be dominated by the transverse con- 
ductivity except in a very narrow range of directions of 
q where the synthetic electric field is almost purely lon- 
gitudinal. 

We are now in a position to combine these results with 
Eq. (|14p and obtain the phonon damping rate in graphene 
over the full range of wavevectors. Whenever the syn- 



8 



thetic electric field has an appreciable transverse compo- 
nent, we find 



(7^< 1 



(49) 



Here, we defined the function f{9) which is equal to unity 
deep in the diffusive regime, and crosses over to f{9) = 
cos^ 30 for transverse phonons and f{6) ~ sin'^ 36* for 
longitudinal phonons'*^. Thus, we find that the phonon 
modes are underdamped at long wavelengths <C 1, but 
become marginal for 3> 1. Inserting numbers, we find 
that Fq w 10~^ajq in this marginal regime so that the 
phonon mode remains well defined. 

Whenever f(9) is close to zero, the damping is domi- 
nated by the longitudinal conductivity. Specifically, this 
happens when q points in the zigzag (armchair) direction 
for longitudinal (transverse) phonons. In these cases, the 
damping exhibits an intermediate q-independent regime 
in between the quadratic and the linear wavevector de- 
pendence. Specifically, we find 



qli/2 
{c/vfY 



9> 7 



(50) 



for very weak intervalley scattering qoiy 3> 1 and 



q'^e/2 



£/2il 



{clvpfq 



9«i 

(51) 



for stronger intervalley scattering qaiy ^ 1- 

Our results for the acoustic (in-plane) phonons of 
graphene are summarized in Figs. [3] and H) For trans- 
verse synthetic fields, the damping rate crosses over from 
a quadratic to a linear dependence on q when q ^ l/£. 
In contrast, when the damping is dominated by longi- 
tudinal fields, there is an intermediate constant regime 
for max{qo,l/^y} < g < max{l, l/((7o^y)2}(l/£). In 
magnitude, the damping rate is the same for longitudi- 
nal and transverse electric fields for small q where we find 
a quadratic dependence on q in both cases. For larger q, 
the damping is much stronger for transverse fields. Deep 
in the ballistic regime q£ ^ 1, where one finds a linear 
dependence on q in both cases, the ratio saturates at ap- 
proximately {vf/cY ~ 10**. 



2. Flexural modes 

We now turn to flexural modes of suspended graphene 
samples which are characterized by out-of-plane fluctua- 
tions /i(r, t). A qualitative difference arises since, by sym- 
metry, h{r,t) appears quadratically in the strain tensor. 
As a result, the coupling between electrons and flexural 
modes is quadratic rather than linear. This implies that 
the dominant damping mechanism involves both the cre- 
ation of a particle-hole pair and a lower-energy flexural 
phonon. 



Gph(q .«) 




Gph(q .«) 



CT(q,00) 



FIG. 5: (Color online). Sketch of the Feynman diagram for 
the damping rate of flexural phonons. 



The corresponding rate can be obtained from Fermi's 
golden rule (or alternatively, the Feynman diagram in 
Fig. [5]). Quantizing the out-of-plane displacements, we 
flnd at zero temperature 



X E EEk^i 

efi>EF a<EF q' 

X 6{ef3 - ea ~ LUc^ + uj. 



VpCr ■ Aq,-c 



oi(q-q')T| 



(52) 



in terms of 



A 



qixq2y + qiyq2x 

qixqix — <lly<l2y 



nn 



2poLJqi y 2pQLJq^ 



.(53) 

Following the approach of this paper, we rewrite this ex- 
pression in terms of the dissipative conductivity Eq. (jl3p . 



Tq = 



1 



2e^hn V 2a 



2po 



X a 



Qxq'y + qyq'x \ f Qxq'y 
qxq'x-qyi'y J , \ qxq'x 



qyqx 
qyq'y 



(54) 



As we have seen above, the conductivity tensor is domi- 
nated by the transverse conductivity so that 



(q,w) ~ CTx(q,Lj) ( % - ^ 



(55) 



At flnite temperature T ^ Wq , the integrand in Eq. (I54p 
is modified to include an extra factor T/fkUqr which is 
a consequence of the existence of a phonon in the final 
state. 

We first focus on ideal graphene membranes for which 
the dispersion of flexural phonons is quadratic, Wq — aq^ 
(where a = y^n/po) . In this case, we find for the angular 
average of the damping rate 



-q 



ql-^l 



,i„2„" q 



(56) 
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FIG. 6: (Color online). Inverse quality factor as function 
of temperature for long wavelength flexural vibrations in 
graphene. The wavevector is |q|^^ = l^m, the mean free 
path is / = lOOnm, and the carrier density is lO^^cm"^. Left: 
No external tension. Right: External tension 7 = 0.02. 



in the low-temperature limit T ^ haq^ and 







(T_ 



3/2 



(57) 



in the high- temperature regime T 3> haq^. Because of 
the low frequency of pure flexural phonons, the vibrations 
can become overdamped at finite temperature. 

If the graphene membrane is under tension 7 (inducing 
a term ~ (V/i)^ in the elastic Lagrangian), the disper- 
sion of the flexural phonons becomes linear at long wave- 
lengths, Wq = Cfq. In this case, we find for the angular 
average of the damping rate 



g£< 1 
1 



(58) 



in the low-temperature regime T <C hcjq and 



P^o" 



T < hcf/l 



(59) 



in the high-temperature T ^ hcfq regime. '^^ 

Actual suspended samples are typically expected to be 
under some degree of tension. The magnitude of phonon 
damping will then depend critically on the sample- 
specific degree of tension, which enters through the mode 
velocity cj — c^/y. It is interesting to point out that un- 
like in-plane modes, the damping of flexural modes is 
temperature dependent in the degenerate limit Ep ^ T, 
even in the absence of electron-electron interactions. 



IV. COULOMB INTERACTION EFFECTS ON 
PHONON DAMPING 

A. Damping by valley Coulomb drag in graphene 

Electron-electron interactions lead to interesting tem- 
perature dependence of the phonon damping which per- 
sists even in clean samples. In fact, it is a well-known 



FIG. 7: (Color online). Temperature dependence of the in- 
verse quality factor Fq/ojq for acoustic phonons in graphene, 
induced by valley Coulomb drag, when damping is dominated 
by longitudinal (left) and transverse (right) synthetic electric 
field. Calculations have been done for an elastic mean free 
path £ = lOOnm, tv ~ 200r, |q| — (l/^m)~^, and a car- 
rier density n = lO^^cm"' 
transverse synthetic fields 



Note the enhancement of Fq for 



effect of spintronics that spin currents decay due to 
electron-electron interactions even in the absence of dis- 
order since unlike charge currents, they are not protected 
by momentum conservation. By analogy, the valley cur- 
rents driven by the electric fields associated with the vec- 
tor potential will dissipate by intervalley Coulomb scat- 
tering. A theory of this effect must go beyond exisiting 
works on Coulomb drag in that one has to account for 
the ac nature of the driving valley-odd electric field. 

This is easily accomplished in the diffusive limit where 
we can amend the hydrodynamic equations to include 
the intervalley electron-electron scattering. While the 
continuity equations remain unchanged, Ohm's laws take 
the form 



1 


1 






1 


- H 


1 






1 — 


r 


TV 


TD J 


\TD 


TV 



J2 



1 

TD 



1 

TV 



Jl 



TV 



-El 
1 

TD 



J2 



(60) 



(61) 



Here, \/td ~ T'^ /Ep (possibly up to logarithmic factors 
from 2kp scattering) is the intervalley scattering rate re- 
sponsible for the drag effect. (Coulomb drag in doped 
graphene is expected to have the same features as in 
other 2D electronic systems.—) It is then evident that 
valley Coulomb drag can be accounted for by evaluating 
the dc conductivity Udc and the diffusion constant D with 
the effective scattering rate 



1 

Toff 



2 

TV 



2 

td' 



(62) 



Thus, the dissipation is dominated by the drag effect (dis- 
order scattering) when T > T* {T <T*), where 



(63) 
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Here, we used that for realistic samples, disorder pre- 
dominantly causes intravalley scattering. For realistic 
parameters in suspended graphene, we estimate T* w 
10 - lOOK. 

The inclusion of valley Coulomb drag allows us to dis- 
cuss the temperature dependence of the intrinsic phonon 
damping in clean (but doped) graphene. The interval- 
ley electron-electron scattering leads to a temperature- 
dependent mean free path ^ HvpEp /T"^ ■ Thus, when 
holding the wavevector q fixed, we find that the damping 
will follow the behavior of the ballistic (diffusive) regime 
for T<T*{T> T*), where 

T* ~ ^EF{hvFq). (64) 

As a result, we find for fixed q and when the damp- 
ing is dominated by the transverse synthetic electric field 
that Fq is constant for temperatures smaller than T* 
(determined by the transverse conductivity in the bal- 
listic regime) and decreases monotonically as oc 1/T^ as 
the system enters the diffusive regime for T > T*. This 
remarkable temperature dependence is a consequence of 
the fact that the dissipation is proportional to the scatter- 
ing time and thus inversely proportional to the scattering 
rate. 

In the presence of disorder, the temperature depen- 
dence also crosses over from constant to monotonically 
decreasing as oc This crossover occurs at T* when 

q£'^ 1 and at T* when g£ <C 1 . When the synthetic elec- 
tric field is purely longitudinal, there can be an interme- 
diate regime (for T > T* in the clean limit) in which the 
damping rate increases as before it crosses over into 
the oc 1/T^ behavior beyond a temperature {vp/cy^^T^- 

B. Damping by valley Coulomb drag in carbon 
nanotubes 

Valley Coulomb drag also affects the damping of 
phonon modes of carbon nanotubes. For the radial 
breathing mode of metallic carbon nanotubes, the effect 
may be significant since phonon damping is dominated 
by the decay of valley currents. Moreover, it is natural 
to expect that Coulomb drag is particularly effective in 
a one-dimensional setting. In contrast, we expect that 
valley Coulomb drag is less significant for the longitudi- 
nal stretching mode where damping is dominated by the 
effects of the scalar potential. In this case, damping by 
decay of valley currents is only a subleading contribu- 
tion. Nevertheless, it is worthwhile to remark that val- 
ley Coulomb drag provides the dominant damping mech- 
anism based on the electron-phonon interaction in the 
strictly clean limit. 

A detailed theory of valley Coulomb drag in carbon 
nanotubes must account for possible electronic correla- 
tion effects associated with the Luttinger liquid nature of 
the electron system. Such a theory is beyond the scope 
of the present paper. 



V. CONCLUSIONS 

We have considered the damping of low-energy 
phonons in carbon nanotubes and graphene originating 
from the electron-phonon interaction. For most phonon 
modes, this damping is closely related to a synthetic elec- 
tric field associated with a strain-induced vector poten- 
tial in the Dirac equation for the electronic properties 
of graphene. We find that it is very instructive to ana- 
lyze phonon damping in terms of these synthetic electric 
fields: (i) Within this approach, phonon damping is a 
direct consequence of Joule heating, (ii) This establishes 
a close relation between phonon damping and the dy- 
namic conductivity which we exploit to derive damping 
rates in the presence of disorder and electron-electron in- 
teractions, (iii) We find rich physics emerging from the 
fact that the synthetic electric field has opposite signs 
in the two valleys. Most prominently, we identify valley 
Coulomb drag as an important dissipation channel which 
leads to unconventional temperature dependence of the 
damping rate. 

Throughout this paper, we have considered idealized 
samples in the sense that we ignored finite size effects and 
electrodes. Clearly, when suspended carbon nanotubes or 
graphene membranes are coupled to electrodes, there will 
be (additional) electronic dissipation taking place in the 
leads even if the nanotube is otherwise perfectly ballis- 
tic. While the physics of this damping is certainly highly 
non-universal, a rough estimate may be obtained from 
our expressions for diffusive electronic dynamics by set- 
ting the elastic mean free path equal to the length of the 
carbon nanotube or the linear dimension of the graphene 
membrane. It is also worthwhile to point out that the 
electron-phonon coupling is not expected to change sig- 
nificantly in multilayer graphene samples. Thus, our re- 
sults should also be applicable in these systems. 

Our results should be of direct relevance to the in- 
tense ongoing experimental efforts to build and explore 
nanomechanical as well as nanoelectromechanical de- 
vices based on graphene nanostructures. We expect the 
electron-phonon interaction to be the dominant source of 
phonon damping whenever the system exhibits a metal- 
lic conductivity. In such systems, our results should be 
valuable by providing upper bounds on the quality factor 
as well as by guiding optimization strategies. 
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Here, we asumed hnear response and Up^^ denotes the 
Fermi-Dirac distribution. Introducing the Fermi surface 
deformation 6u((j)) by Aup — S{ep — fj,)5iy{(j)) and intro- 
ducing the angle = Z(q, p) = ^(E, p), we find 



APPENDIX A: DYNAMIC CONDUCTIVITY OF 
GRAPHENE IN THE BALLISTIC REGIME 

Here, we sketch the derivation of the longitudinal con- 
ductivity in the ballistic regime, including the effect of 
intervalley disorder scattering. We start from Boltzmann 
equations with valley electric field E for the two valleys, 



, rJ.^ eEvF COS (f) - :^{Siy{(t>)) 

H^) = —. \ , 1 ■ (A4) 

—loj + ivpq cos (p + — 



In the absence of interlayer scattering, we find for the 
density p = ev J {d(j)/2'K)Sh'{(l>) the conventional result 



eE 



{-ioj + ivp ■ q + — ■ Vp)nW 

= fi^E^(^p-V)(-?^-«[>^») (Al) 



eE 



-luj + iVp • q 



2 Vp)<' 



^X:^(^p-v)(-?-ni,^») (A2) 



where rip denotes the distribution function in valley j. 
Here we include the intervalley disorder scattering which 
is required to obtain results which match the diffusive re- 
sults. Taking the difference between these two equations, 
we obtain an equation for the odd distribution function 



Aup ~ n 



(1) 



,(2)) 



{-iuj + ivp ■ q H )Anp + eE • VpTi^'' 



(cq) 



(A3) 



p' 



p=—E. 

iq 



(A5) 



From the expression for the current, j = 
evvp ^ {d(j)/ 2t:) cos (j)5v{(j))^ we obtain the (dissipative) 
longitudinal conductivity 



9 2 

e^^/w 



(A6) 



Eq. (|A4p shows that in the presence of intervalley scat- 
tering, there is an additional contribution (proportional 
to l/ry) to the current. In this term, we can evaluate 
{Si'{(f>)) using the result for Si'{(j)) ~ eE/iq in the absence 
of intervalley scattering. Noting that we may also replace 
cos (f) in the first term by uj/vpq, we find that the second 
term dominates as long as q ^ [(9o^y)^]^^ while the first 
term dominates for q ^ [{qQf-v)C]~^ ■ Evaluating the lon- 
gitudinal conductivity then gives the results quoted in 
Eq. gZD. 
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We note that the description of the synthetic fields as elec- 
tric and magnetic fields is motivated by their effect on the 
Dirac carriers. At the same time, the dynamics of these 
fields is entirely controlled by elasticity theory which has 
little resemblence to electrodynamics. 
The crossover occurs at q ^ 1/i when qoiv S> 1 and at 
q ~ l/£v when qoi <C 1. A more accurate expression for 
the function f{9) can be readily derived by including both 
the longitudinal and the transverse components of the syn- 
thetic electric field in the calculation of the damping rate. 
In principle, there are also intermediate cases in which the 
dispersion is linear at long wavelengths due to tension but 
becomes quadratic at larger q. Our results can be extended 
to this situation but we refrain from giving explicit results. 



